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摘 要 :有 限 元 分 析 方 法 中 ,采用 三 维 实体 单元 对 焊接 管 结 构 进 行 离散 是 最 为 精确 的 建 模 方法 ,但 
对 三 维 实体 模型 分 析 时 会 划分 大 量 单元 ,导致 计算 时 间 史 长 。 工 程 中 常 把 焊接 管 结构 简化 为 刚 加 
模型 , 即 刚性 连接 的 杆 系 系统 。 这 种 建 模 方式 忽略 了 节点 的 局 部 柔 度 , 过 高 地 估计 了 整体 结构 的 刚 
度 , 使 计算 结果 偏 于 危险 。 为 解决 上 述 问题 ,在 节点 处 引入 虚拟 粱 单元 ,使 其 等 效 模拟 节点 的 局 部 
变形 ,还 原 支 管 真实 长 度 。 为 了 验证 节点 柔 度 对 焊接 管 结构 的 影响 , 对 典型 的 焊接 管 结构 缩 尺 , 进 
行 了 推荐 试验 。 分 别 采 用 沉 单 元 模型 \ 刚 架 模型 及 引入 虚拟 梁 单 元 的 简化 模型 进行 有 限 元 模拟 ,并 
-二 与 试验 结果 对 比分 析 。 结 果 表 明 , 刚 架 模型 会 过 高 的 估计 焊接 管 结构 的 刚度 ,而 简化 模型 计算 结果 
) 与 过 模型 .试验 结果 较为 吻合 。 
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Simplified nonlinear pushover model of welded tubular structure 
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Abstract: Among the finite element (FE) analysis methods, using three-dimensional solid elements to 
discrete the welded tubular structure is the most accurate modeling method. The 3D solid model involves a 
“huge number of units during the analysis process ,which results in time-consuming calculation. Therefore ， 
an accepted alternative is simplifying the welding tubular structure to a rigid frame model in engineering, 
that is ,a rigidly connected rod system. However ,the local joint flexibility ( LJF) is neglected and the rigid- 
ity of the structure is overestimated in the simplified steel frame model,which is likely to lead to slightly 
dangerous calculation results. A fictitious beam element (FBE) is introduced into the joint to simulate the 
local deformation of the joint and restore the true branch tubular length ,which can solve the above prob- 
lems effectively. A push-over testing on the scaled typical welded tubular structure was carried out to verify 
the influence of joint flexibility on welded tubular structure. The shell element model ,the rigid frame mod- 
el,and the mentioned simplified model were adopted for FE simulation ,respectively ,and the results were 


compared with FE and experimental results. The results show that the rigid frame modeling will overesti- 
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mate the stiffness of the welded tubular structure, while the simplified modeling results meet well with 


those of the shell modeling and testing. 
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采用 相 贯 节点 的 焊接 管 结构 有 着 造型 简洁 、 质 
量 较 轻 .连接 方便 等 特点 ,因而 广泛 应 用 于 石油 工 
业 ,城市 建筑 .电力 工程 等 方面 。 对 复杂 焊接 管 结构 
进行 力学 分 析 时 常 将 其 简化 为 刚 架 模型 。 因 为 主管 
管 壁 较 薄 ,主管 的 径 向 刚度 远 小 于 支管 的 轴 向 刚度 ， 
所 以 节点 处 主管 壁 更 易 受 载 变形 ,这 种 现象 称 为 局 
部 柔 度 。 刚 架 模型 忽略 了 节点 的 局 部 柔 度 ,使 整体 结 
构 的 刚度 偏 高 "。 针 对 这 一 问题 ,以 刚 架 模型 为 基 
础 ,在 节点 处 引入 虚拟 梁 单 元 以 模拟 节点 局 部 变形 。 

Boukamp 等 中 最早 研究 节点 局 部 柔 度 对 海洋 平 
侣 的 影响 ,对 用 于 建立 节点 柔 度 模型 的 分 析 技 术 进 
征 醒 研究 ,并 将 这 些 技术 应 用 于 平台 整体 中 。 结 
表明 节点 柔 度 对 整体 结构 有 着 显著 的 影响 , 且 在 构 
件 长 度 短 于 公称 直径 的 情况 下 影响 最 为 明显 。Mir- 
tdhehi 等 中 建立 了 两 个 平台 的 有 限 元 模型 ,研究 节点 
革履 对 固有 频率 的 影响 ,进行 了 非 线性 静 力 与 动力 
分 击 。 结 果 表明 节点 柔 度 对 海上 结构 的 非 线性 动力 
特 钥 有 着 显著 影响 。 为 了 系统 地 研究 钢管 参数 对 节 
点 局 部 柔 度 的 影响 ,国内 外 众多 专家 学 者 针对 各 种 
类 型 的 节点 进行 模型 试验 、 半 解析 法 [ss 有 限 元 
模 筷 ”中 得 到 了 大 量 的 节点 局 部 柔 度 参 数 公 式 。 文 


限 元 仿真 ,通过 试验 验证 了 模型 分 析 方 法 的 准确 性 ， 
对 有 限 元 结果 进行 回归 分 析 , 得 到 了 弯 矩 /转角 参数 
公式 。 该 公式 可 以 很 好 地 描述 钢 框 架 节 点 受 弯 产生 
的 非 线 性 变形 。 局 部 柔 度 参数 公式 不 仅 有 利于 研究 
其 对 焊接 管 整 体 结构 的 力学 影响 ,而 且 使 节点 局 部 
柔 度 等 效 单元 的 实际 应 用 更 加 方便 。 

朱 静 "在 主管 和 任意 二 支管 间 建 立 虚拟 单元 ， 
以 等 效 节 点 局 部 柔 度 ,并 举例 对 比 该 模型 与 传统 刚 
架 模型 的 力学 性 能 。 结 果 表 明 两 种 模型 力学 性 能 
着 较 大 差异 。Golafshani 等 基于 柔性 矩阵 和 
Fessler 经 验 公式 趾 提 出 了 一 种 新 的 局 部 柔 度 单 元 ， 
使 用 两 个 连接 节点 将 主管 壁 变形 合并 到 整体 刚度 所 
阵 , 但 模型 的 应 用 较为 复杂 且 仅 限于 弹性 范围 内 。 
安利 强 等 “将 相 贯 节点 局 部 柔 度 公式 应 用 于 插 板 
节点 ,把 连接 螺栓 超过 两 个 的 节点 处 理 为 半 刚 性 节 
点 ,分 析 了 刚度 因子 取 值 对 杆 件 轴 力 与 节点 位 移 的 
影响 ,结果 表明 刚度 因子 取 值 对 二 者 影响 显著 。 李 
东升 ”在 支管 端 部 附加 弹 自 和 扭 得 ,使 其 蔡 代 节 点 
变形 。 使 用 虚 功 原理 推导 出 了 考虑 节点 柔 度 梁 单元 
的 刚度 矩阵。 针对 海上 平台 典型 的 平面 焊接 管 结构 
利用 两 种 建 模 方式 进行 推 履 分 析 , 结果 表明 考虑 节 


献 癌 ] 进行 了 载荷 作用 下 T、Y 塑料 节点 的 变形 实 
验 5 利 用 实验 结果 得 出 节点 在 轴 力 、 弯 和 矩 作用 下 局 部 
骤 厦 的 参数 公式 。 陈 伯 真 等 "根据 梁 的 弯曲 理论 得 
到 QDY 节点 局 部 柔 度 公 式 , 计 算 了 8 组 节点 的 局 部 
柔 度 , 计 算 结 果 与 Fessler 参数 公式 计算 结果 较为 
符合 。 文献 [9] 利 用 有 限 元 方法 对 节点 单 参 分 析 ， 
将 计算 结果 进行 多 元 回归 分 析 得 到 了 TY、K 节点 
的 参数 公式 。 文 献 [11] 在 两 种 轴 向 载荷 作用 下 对 
81 个 PE 模型 进行 了 模拟 ,根据 模拟 结果 比较 空间 
DK 节点 与 K 节点 的 局 部 柔 度 因子 (fi1s)。 结 果 表 
明 空间 节点 与 平面 节点 的 局 部 柔 度 差异 很 大 ,针对 
这 一 问题 开发 了 两 个 单独 的 参数 公式 来 确定 受 两 种 
轴 向 载荷 的 空间 DK 节点 中 的 fijs 值 。 陈 映 等 中 建 
立 了 63 个 型 弯 管 锥 形 节点 有 限 元 模型 ,研究 了 径 
厚 比 、 直 径 比 、 主 管 与 支管 夹 角 等 因素 对 节点 轴 向 
刚度 的 影响 ,并 通过 拟 合 得 到 弹 塑性 范围 内 的 节点 
轴 向 刚度 因子 与 轴 力 因子 的 函数 关系 。 李 明 吴 
等 (中 针对 几 种 典型 钢 框 架 节 点 的 抗 弯 性 能 进行 有 


点 柔 度 的 梁 单元 程序 具有 良好 的 适用 性 。 王 骼 
在 节点 处 引入 虚拟 梁 单 元 以 模拟 节点 的 局 部 变形 ,对 
比分 析 了 毅力 作用 下 平面 TY 节点 管 机 架 在 弹性 和 
弹 塑 性 范围 内 利用 3 种 建 模 方 式 所 得 到 的 计算 结果 ， 
并 进行 了 参数 分 析 。 结 果 表 明 考 虑 节点 柔 度 的 简化 
模型 可 以 很 好 地 模拟 TY 节点 管 桥架 的 节点 柔 度 。 

已 有 的 考虑 节点 柔 度 的 简化 分 析 方 法 大 都 限制 
在 弹性 范围 内 ,王皓 "” 提出 的 简化 模型 通过 与 壳 单 
元 模型 对 比 ,验证 了 该 方法 在 平面 T、Y 节点 管 析 架 
问题 中 的 适用 性 。 本 研究 在 虚拟 梁 单 元 模型 基础 
上 ,针对 空间 焊接 管 结构 进行 推荐 分 析 , 以 研究 虚拟 
梁 单 元 模型 在 分 析 静 力 非 线性 问题 时 的 可 行 性 。 


1 节点 局 部 柔 度 的 定义 和 虚拟 梁 单 元 
模型 


1.1 节点 局 部 柔 度 定义 
焊接 管 结构 支管 受 载荷 作用 ,节点 处 的 主管 管 
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壁 会 产生 局 部 变形 ,那么 节点 处 便 具 有 了 局 部 柔性 。 各 


入 与 主管 轴线 垂直 的 虚拟 粱 单元 ,如 图 2 所 示 。 
对 于 任意 节点 ,节点 局 部 柔 度 定 义 为 ‖ 上 可 六 站 元 
A FE=-=- 二 -= 一 -一 -一 -= [7 
AX 二 万 (1) 
Oo (Ga) 普通 T 节 点 的 简化 
Olps = 37 (2) D2 
. | 0 
5 (3) ED 
OPB ~ M, 


式 中 :4 为 轴 力 作用 下 节点 处 主管 管 壁 因 变形 产生 
的 垂直 于 主管 轴线 方向 的 位 移 ;P 为 轴 力 ;p, 为 平面 
内 弯 矩 作用 下 主管 管 壁 的 面 内 转角 ;M, 为 主管 与 支 
管 所 在 平面 内 的 弯 算 ;po 为 平面 外 弯 矩 作用 下 主管 
管 壁 的 面 外 转角 ; Wo 为 主管 与 支管 所 在 平面 外 变 
和 矩 。 如 图 1 所 示 , 在 主 .支管 相 贯 线 上 ,a .2 为 冠 点 ， 
c 可 为 鞍点 ;e 为 主 \ 支 管 轴线 交点 ,人 g 为 过 。 点 且 
重重 于 主 支管 轴线 所 形成 平面 的 直线 与 主管 外 辟 
的 两 个 交点 。A, .A; .A。 .As、Aj、As 分 别 表 示 对 应 点 
玲 放 于 主管 轴线 方向 的 位 移 ,其 中 人 g 点 位 移 由 主 
管 酌 弯曲 引起 ,ac sd 点 的 位 移 由 主管 弯曲 和 节点 
处 者 管 壁 变形 引起 ,而 主管 壁 的 变形 可 以 由 各 种 载 
荷 习 起 。 


图 1 YY 市 点 示意 医 
Fig.1 View of a tubular Y-joint 


.二 负 力 作用 下 节点 处 主管 管 壁 垂直 于 主管 轴线 方 


向 前 位 移 定义 为 
| a Aj+A, 4 
和。 4 D> 9 
平面 内 弯 矩 作用 下 主管 管 壁 的 面 内 转角 定义 为 
py = arotan en (5) 


式 中 :0 为 主管 与 支管 的 夹 角 ;du 表示 文 管 外 径 。 同 
理 ,平面 外 弯 矩 作用 下 主管 管 壁 的 面 外 转角 定义 为 
| A, -A |sing 


po = arctan 


1.2 虚拟 梁 单 元 模型 

传统 刚 架 模型 将 焊接 管 简 化 为 其 轴线 处 的 梁 单 
元 , 梁 单 元 间 刚 性 连接 。 因 主管 管 壁 厚度 远 小 于 文 
管 长 度 , 则 当 支 管 端 部 在 轴 力 作用 下 时 ,主管 径 向 刚 
度 便 远 小 于 支管 轴 向 刚度 ,节点 处 主管 壁 的 变形 便 
几乎 垂直 于 主管 轴线 ” , 则 可 在 刚 架 模 型 节点 处 引 


(b) 普通 Y 节 点 的 简化 

妈 2 简化 模型 的 简化 过 程 
Fig.2 Simplified process of simplified model 

虚拟 梁 单 元 的 柔 度 有 其 材料 属性 和 几何 属性 共 
同 确定 。 在 柔 度 确 定 的 情况 下 ,虚拟 梁 单 元 的 弹性 
模 量 k 可 以 选取 任意 值 ,因为 可 以 通过 调整 几何 参 
数 来 符合 确定 的 柔 度 值 。 实 际 模型 为 闭口 截面 构 
件 , 抗 扭 刚度 较 大 ,扭转 变形 较 小 ,所 以 虚拟 梁 单 元 
不 考虑 平面 外 弯 矩 作用 。 虚 拟 梁 单元 可 以 等 效 为 固 
定 在 主管 上 的 基 臂 梁 ,依据 材料 力学 知识 及 节点 和 柔 
度 的 定义 可 以 得 到 虚拟 梁 单 元 的 截面 积 4 和 截面 惯 
性 矩 7。 


/三 Wy 
” 2 有 .pm 


虚拟 梁 单元 截面 在 满足 4 和 7 的 条 件 下 可 以 取 
任意 形状 。 为 保持 整体 结构 的 一 致 性 ,虚拟 梁 单 元 
假定 为 圆 钢 管 结构 ,其 截面 积 4 和 截面 惯性 矩 1 为 
Tdi —-d) T(d -di) 

4 2 4 
由 式 (6) 可 得 虚拟 染 单 元 外 径 d 和 内 径 心 , 即 


2 焊接 管 结构 静 力 推 履 试验 


2.1 试 件 设计 
焊接 管 结构 试 件 包含 了 节点 下 节点 、 空 间 TY 
节点 这 3 种 常见 的 节点 类 型 ,结构 的 设计 符合 
《GB 50017 一 2017 钢 结构 设计 标准 规范 》。 试 件 的 
几何 结构 如 图 3 所 示 。 主 、 文 管 尺寸 如 表 1 所 示 。 
表 1 主 、 支 管 尺寸 


Tab.1 Dimensions of chords and braces 


A= 


名 称 主管 横 支 管 斜 支管 
外 径 /mm 114 60 33.5 
厚度 /mm 3 2 1.5 


7 
T 
1 
1 
1 
1 
1 
1 
1 
1 
1 


和 


~ 


(a) 推 覆 试 件 正 视 医 (b) 推 履 试 件 俯视 
图 3 推荐 试 件 的 几何 斥 才 
Fig.3 Geometry of test specimen 
选用 Q235A 作为 试 件 材 料 。 分 别 截取 焊接 管 
架 顶 板 、 主 管 . 横 支 管 \, 斜 支管 材料 制作 试 件 且 拉 伸 
试 件 与 焊接 管 架 是 来 自 同 一 批 次 的 钢材 , 借 此 消除 
厚度 、 曲 率 对 材料 属性 的 影响 。 严 格 按照 
《GB/AT 228. 1 一 2010 金属 材料 的 室温 拉 伸 试验 》 设 


绩 硅 均值 以 减 小 偶然 误差 。 各 组 试 件 的 材料 属性 数 


据 如 表 2 所 示 。 
OO 表 2 试 件 的 材料 属性 
CD Tab.2 Material properties of test specimen 


材料 属性 顶板 主管 横 支 管 。 斜 支管 


模 量 /GPa 218 220 231 235 
服 强度 /MPa 238 245 240 220 
强度 /MPa 391 385 369 363 
率 /% 25 21 24 22 


2, 公 试验 及 测量 装置 

< 依靠 一 个 1000 kN 水 平 作 动 器 为 试验 施加 推荐 
力 回 作 动 器 一 端 固定 至 反 力 墙 , 另 一 端 与 连接 件 形 
焊接 管 斌 
人 


脚 焊接 在 框 型 底板 上 ,底板 与 底座 通过 螺栓 连 
接 实 现 固 定 目 的 ,加 载 装置 示意 图 如 图 4 所 示 。 推 
窗 试 件 为 对 称 结构 ,试验 加 载 尽量 接近 中 线 , 两 榴 框 
架 在 加 载 过 程 中 扭转 的 作用 很 小 ,可 以 忽略 。 从 上 
至 下 设置 了 5 个 高 精度 位 移 计 , 依 次 测量 作 动 咒 、 顶 
板 、 横 支管 节点、 主管 的 位 移 ,如 图 5 所 示 。 


图 4 推荐 试验 加 载 示 意图 


Fig.4 Loading scheme of pushover test 
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图 5 ”位移 计 安装 位 置 
Fig.5 Placement of LVDTs 

2.3 试验 过 程 

以 控制 加 载 位 移 的 方式 对 推 覆 试 件 进行 加 载 ， 
直至 观察 到 明显 的 变形 或 破坏 便 停 止 加 载 。 试 验 中 
与 作 动 器 平行 的 2 个 K 节点 均 发 生 了 严重 的 破坏 ， 
一 处 节点 斜 支管 发 生 明显 的 弯曲 变形 ,主管 向 下 
弯曲 , 另 一 处 K 节点 的 主 、 支 管 焊接 处 发 生 断 裂 , 顶 
板 与 主管 弯曲 。 图 6 和 图 7 分 别 给 出 了 试 件 的 整体 
变形 和 局 部 变形 情况 。 


图 6 试 件 的 整体 变形 


Fig.6 Global deformation of test specimen 


局 


节点 焊接 处 断裂 国 国 


了 节点 支管 处 弯 


(@) K 节 点 焊接 处 断裂 (b)K 节 点 主 、 支 管 弯 曲 变形 
图 7 试 件 的 局 部 变形 


Fig.7 Local deformation of specimen 


3 ”虚拟 梁 单元 几何 尺寸 及 材料 属性 


节点 局 部 柔 度 可 以 利用 已 有 参数 公式 得 到 ,也 
可 对 单一 节点 有 限 元 分 析 得 到 。 因 为 空间 TY 节点 
参数 公式 尚未 提出 ,本 着 一 致 性 的 原则 ,对 试验 涉及 到 
的 TT 空间 TY 节点 均 进 行 有 限 元 分 析 得 到 局 部 柔 
度 。 对 各 节点 支管 施加 合适 的 轴 力 和 面 内 弯 矩 载 葵 ， 
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使 其 在 弹性 范围 内 有 着 尽量 大 的 变形 ,以 减 小 误差 。 
3.1 几何 尺寸 

以 TT 节点 为 例 ,建立 TT 节点 三 维 壳 单元 模型 ， 
主管 长 800 mm ,支管 均 为 150 mm , 除 加 载 端 , 主管 及 
支管 端 部 均 设 置 成 固定 端 。 在 支管 端 部 施加 15 kN 
的 轴 力 ,如 图 8 所 示 。 加 载 过 程 始 终 处 在 弹性 阶段 。 
采用 S4R 单元 对 模型 进行 网 格 划 分 。 输 出 a、b、e、 
df 点 位 移 , 结 合式 (2) 得 到 轴 力 作用 下 节点 处 主管 
管 壁 因 变形 产生 的 垂直 于 主管 轴线 方向 的 位 移 ,再 
代入 式 (1) 得 到 TT 节点 轴 向 柔 度 。 利 用 相同 模型 ， 
在 支管 端 部 施加 500 N . m 的 弯 抢 ,其 余 主 支管 端 
面 设 置 成 固定 端 , 如 图 9 所 示 。 加 载 过 程 始 终 处 于 
模型 的 弹性 阶段 。 


图 8 TT 市 点 模型 受 轴 力 作用 


Fig.8 TT-joint under axial compression 
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图 9 TT 节点 模型 受 弯 矩 作用 
Fig.9 TT-joint under in-plane bending 

输出 ab 两 点 位 移 代 和 人 式 (3) ,可 以 得 到 TT 节 
点 面 内 弯 矩 作用 下 的 面 内 转角 ,再 将 面 内 转角 代入 
式 (1) ,得 到 平面 内 论 曲 柔 度 。 取 所 有 虚拟 粱 单元 
的 弹性 模 量 与 主管 一 致 ,为 220 CPa ,将 2 个 柔 度 值 
代入 式 (5), 则 可 以 得 到 虚拟 梁 单 元 的 截面 积 和 截 
面 惯性 矩 。 同 理 , 可 求解 其 他 节点 的 虚拟 梁 单 元 几 
何 尺 寸 。 各 节点 载荷 、 柔 度 及 虚拟 梁 截面 几何 参数 
如 表 3 所 示 。 表 中 :TY-1/2 中 的 1 表示 和 斜 支管 ,2 表 
示 与 主管 垂直 的 文 管 ; 表 中 节点 的 轴 力 和 弯 矩 理论 
可 以 取 弹 性 范围 内 的 任意 一 组 载荷 ,所 得 到 的 柔 度 
值 均 相 同 ,但 为 减 小 误差 ,选择 弹性 范围 内 的 较 大 
载 信 。 


表 3 各 市 点 的 载荷 . 柔 度 及 虚拟 梁 截 面 几何 参数 


Tab.3 Load,flexibility and geometric parameters of the fictitious beam element at joints 


a 载 向 邓 度 虚拟 梁 单 元 截面 
F/ M/ 6aAx/ Slpp/ 
Ga A /mm IT /mm’ di /mm dy/mm 
加 kN (CN mm) (mm . N-!) (rad.N- mm-1) 
T= 15 500 9.09E-06 6.16E-09 28.50 42 060 108.74 108.57 
= 10 100 5.18E-06 1.66E-08 18.62 8 215 59.52 59.32 
Two-planar TY-1 1 150 7.46E-06 1.43E-08 24.56 18 144 76.99 76.78 
Two-planar TY-2 1.2 500 3.36E-06 5.67E-09 77.11 45 695 69.21 68.50 
3.2 材料 属性 


随 着 载荷 的 增加 ,材料 进入 塑性 阶段 ,载荷 与 位 
移 不 再 是 线性 关系 。 节 点 柔 度 值 发 生变 化 ,不 能 
徘 弹性 阶段 节点 参数 公式 求 出 ,可 以 通过 对 节点 的 
有 限 元 分 析 得 到 。 与 求解 虚拟 梁 单 元 几何 尺寸 类 
似 , 以 空间 TY 节点 为 例 。 建 立 节 点 的 三 维 充 单 元 
模型 ,主管 长 度 为 800 mm ,支管 长 度 均 为 150 mm。 
研究 垂直 文 管 受 轴 力 作用 的 空间 TY 节点 时 ,在 一 
个 与 主管 垂直 的 支管 端 部 施加 10 mm 的 轴 向 位 移 载 
集 , 其 余 支 管 及 主管 端 部 均 设置 为 固定 端 。 图 10 为 
空间 TY 节点 的 受 载 后 变形 , 受 载 后 节点 处 变形 明 
显 且 处 于 塑性 阶段 。 


图 10 


空间 TY 节点 的 受 载 后 变形 
Fig. 10 Deformation of two-planar tubular TY-joint 
输出 点 ab、c.d、f 位 置 处 的 位 移 及 加 载 点 的 反 
作用 力 , 结 合式 (2) 可 求 得 轴 力 作用 下 点 处 垂直 
于 主管 轴线 方向 的 柔性 变形 。 应 变 与 位 移 的 关系 为 


Psin0 
0 = 一 一 一 


A = Al/lo 


(10) 


真实 应 力 ox 与 真实 应 变 er 定义 为 


or= -nor=0(1 -6) (11) 
0 


将 载荷 位移 利用 式 (8) 转化 得 到 应 力 、 应 变 。 
将 得 到 的 应 力 、 应 变 利 用 式 (9) 转化 成 虚拟 梁 单元 
的 真实 应 力 o1、 真 实 应 变 e+。 图 11 为 载荷 -位 移 曲 
线 , 图 12 为 虚拟 梁 单 元 的 真 应 力 - 真 应 变 曲线 。 虚 
拟 梁 材料 的 塑性 应 变 es 定义 为 


ep=er -Om/b (12) 
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图 12 ”真实 应 力 -应 变 曲线 
Fig. 12 True stress-strain curve 
得 到 的 真实 应 力 与 塑性 应 变 便 是 虚拟 梁 单 元 材 
料 的 塑性 属性 。 其 他 节点 处 虚拟 梁 单 元 的 本 构 关系 
求解 方法 与 之 同 理 ,不 再 一 一 歼 述 。 


4 推荐 试验 的 有 限 元 模拟 


无 论 使 用 壳 单 元 还 是 实体 单元 进行 有 限 元 分 
析 ,都 需要 对 模型 和 载荷 进行 简化 ;因为 ,完全 建立 
连接 件 、 作 动 器 、 反 力 墙 的 模型 不 仅 增 加 分 析 负 担 ， 
研究 的 意义 也 非常 有 限 。 由 于 项 板 较 薄 , 仅 为 
15 mm , 顶板 的 刚度 小 于 连接 件 。 整 个 推荐 过 程 连 
接 件 均 未 观察 到 宏观 的 变形 ,所 以 可 将 连接 件 视 为 
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刚体 ;连接 件 与 顶板 螺栓 连接 可 靠 , 亦 可 将 与 连接 件 
连接 的 部 分 顶板 视 为 刚体 。 
4.1 壳 单 元 模型 

对 推荐 试 件 建立 壳 单 元 模型 , 材料 属性 已 由 拉 
伸 试 验 得 到 ,各 部 分 管 结构 具体 的 弹 塑性 数值 依次 
输入 材料 属性 模块 中 。 依 据 上 面 的 简化 方法 将 顶板 
中 心 350 mm x350 mm 正方 形 区 域 设置 为 刚体 , 控 
制 点 在 其 形 心 上 方 150 mm 处 。 限 制 4 个 柱 脚 底面 
圆 全 部 的 空间 位 移 和 转角 ,使 其 成 为 固定 端 。 因 为 
推荐 试验 中 位 移 计 1 所 测 得 的 位 移 数值 最 大 为 
50 mm( 横 向 位 移 为 50 mm 时 ,K 市 点 焊接 处 断裂 ， 
停止 试验 ) ,那么 在 控制 点 处 施加 50 mm 的 位 移 载 
荷 。 壳 模型 推 履 试验 如 图 13 所 示 。 可 以 看 出 ,顶板 
和 支管 的 变形 与 试验 基本 吻合 ,这 验证 了 建 模 方式 
的 正确 性 。 


a| 


(a) 过 模型 (b) 推荐 变形 区 
图 13 壳 模 型 推荐 试验 

Fig. 13 Pushover test of shell model 

4.2 刚 架 模型 
同样 利用 壳 单 元 建立 顶板 ,焊接 管 由 轴线 处 的 

梁 单 元 建立 ,将 两 者 绑 定 在 一 起 ,得 到 刚 架 模型 。 刚 
体 及 控制 点 设置 与 过 单元 模型 相同 , 梁 单 元 截面 几 
何 尺 寸 即 试 件 的 实际 横 截 尺寸 , 刚 架 模 型 的 边界 条 
件 和 载荷 也 与 壳 单 元 模型 相同 。 刚 架 模 型 推 覆 试验 
如 图 14 所 示 。 


(b) 推荐 变形 图 
图 14 刚 架 模型 推荐 试验 


Fig. 14 Pushover test of rigid frame model 


(a) 刚 架 模型 
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4.3 引入 虚拟 梁 单元 的 简化 模型 


一 
nw (Orfskiss+sn 办 一 一 


架 模型 相同 。 虚 拟 粱 单元 的 材料 属性 及 截面 几何 尺 上 | | mse 
十 前 文 已 经 给 出 ,将 其 直接 引入 属性 模块 即 可 。 未 ,| ZN 
放 缩 的 简化 模型 变形 如 图 15(b) 所 示 。3 种 模型 的 | | = T 
网 格 信息 及 分 析 效 率 如 表 4 所 示 。 因 为 3 种 模型 的 | 
顶板 均 采 用 种 子 尺寸 为 20 的 S4R 单元 ,所 以 表格 中 四 入 化 们 于 (0) 容 变 形 攻 
统计 的 网 格 信息 不 考虑 顶板 。 0 


Fig.15 Pushover test of simplified model 
表 4 3 种 模型 的 网 格 及 分 析 时 间 对 比 


Tab.4 Comparison of FE mesh and computing time for three models 


模型 类 别 种 子 尺寸 局 部 加 密 单元 类 型 


站 节点 数目 单元 总 数 分 析 时 间 / s 
壳 模 型 20 15 S4R 16 013 16 242 611 
于 = 刚 架 模型 100 100 B33 184 200 155 
全 简化 模型 100 100 B33 192 203 102 
全 线 。 图 16 为 推荐 试验 和 3 种 有 限 元 分 析 方式 在 4 
ss 结果 分 析 


个 位 移 计 安 装 位 置 的 载荷 -位 移 曲 线 。 位 移 计 2 位 
于 顶板 位 置 , 在 推 覆 过 程 中 因 顶 板 变 形 而 滑脱 , 且 经 
= 输出 3 种 有 限 元 模型 位 移 计 安装 位 置 的 位 移 以 。” 过 模拟 1、2 位 移 计 处 的 载荷 位 移 曲 线 结果 几乎 重 
制 点 处 的 反作用 力 ,得 到 相应 的 载荷 -位 移 曲 。 合 , 所 以 不 考虑 位 移 计 2 处 的 结 
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16 ” 试 件 不 同位 置 的 载荷 -位 移 曲 线 


Fig. 16 Load-displacement curves at different locations 


由 图 16 可 知 , 有 限 元 计算 结果 与 原始 试验 结 间 和 


相 比 ,极限 载荷 吻合 较 好 ,但 刚度 误差 较 大 。 其 原因 
为 连接 匀 与 顶板 的 螺栓 为 M36 高 强 螺栓 ,螺栓 孔 直 
径 39 mm。 因 此 ,螺栓 连接 受 剪 切 作 用 可 能 在 板 间 
产生 最 大 为 3 mm 的 相对 滑动 。 分 析 试 验 数 据 发 现 : 
位 移 计 1 位 于 作 动 涡 位 置 ,经 有 限 元 计算 可 知 ,顶部 
连接 匀 处 所 受 剪 力 较 大 ,会 产生 较 大 相对 滑动 ;底部 
柱 脚 所 受 剪 力 较 小 ,产生 滑动 较 小 ; 测 得 的 刚度 受 项 
部 滑动 和 底部 柱 脚 处 滑动 的 共同 影响 ,刚度 低 于 有 
限 元 结果 。 位 移 计 3、4 和 5 处 测 得 的 刚度 只 会 受到 


一 刚 架 模型 


一 0 一 试验 
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位 移 /mm 
(a) 位 移 计 1 
100 


底部 柱 脚 处 滑动 的 影响 ,对 刚度 的 影响 较 小 。 

现 将 螺栓 孔 直 径 与 螺纹 大 径 的 差 值 作为 最 大 滑 
动 位 移 剔 除 , 剔 除 滑 动 位 移 后 的 载荷 -位 移 曲线 如 图 
17 所 示 。 可 以 看 出 , 刚 架 模型 在 相同 位 移 处 的 载 葵 
自始至终 明显 大 于 其 他 两 种 建 模 方式 及 试验 的 载 
集 。 简 化 模型 与 序 模 型 在 弹性 阶段 吻合 较 好 ,塑性 
阶段 载荷 略 低 于 充 模 型 ,分析 原 因应 该 是 虚拟 染 单 
元 的 弯曲 变形 会 使 节点 产生 沿 主管 方向 的 变形 ,使 
整体 柔 度 增 加 。 无 论 是 简化 模型 还 是 过 模型 都 与 试 
验 结果 相符 合 。 
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图 17 数据 处 理 后 的 载荷 -位 移 曲线 


Fig. 17 Load-displacement curves after processing data 
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针对 工程 中 常用 的 刚 架 模 型 会 忽略 节点 局 部 柔 
度 , 导 致 结构 整体 柔 度 降 低 这 一 问题 ,在 虚拟 梁 单 元 
模型 基础 上 将 模型 应 用 于 包含 空间 节点 的 空间 问 
题 。 通 过 对 节点 有 限 元 分 析 的 方式 获取 虚拟 染 单 元 
的 截面 尺寸 及 材料 属性 。 对 比 3 种 有 限 元 方法 与 试 
验 获 得 的 载荷 位 移 曲 线 可 以 得 出 以 下 结论 。 

1) 简 化 模型 在 弹性 变形 阶段 与 过 模型 一 致 性 较 
好 ,塑性 阶段 刚度 略 低 于 壳 模 型 ,两 种 建 模 方 法 得 到 


的 结果 均 与 试验 结果 相符 合 。 这 说 明 简化 模型 能 够 
较 好 地 模拟 焊接 管 结构 的 弹 塑性 推荐 ,具有 在 实际 
工程 中 运用 的 价值 。 

2) 简化 模型 相 较 于 壳 模 型 可 以 大 大 提高 计算 效 
率 。 对 推荐 试验 的 有 限 元 分 析 可 知 , 使 结果 收敛 所 
需 的 节点 数 为 元 模型 的 1.20% ,单元 数 为 这 模型 的 
1.25% ,计算 时 间 为 元 模型 的 16.69% 。 因 简化 模型 
顶板 的 网 格 类 型 和 尺寸 与 却 模 型 相同 ,顶板 的 网 格 
数 为 简化 梁 单 元 的 19.45 倍 ,所 以 , 当 模型 较为 复杂 
或 板结 构 能 够 简化 成 刚体 时 ,简化 模型 的 计算 效率 
能 够 得 到 更 大 提升 。 
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